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Abstract- The present paper is aimed at to study the reflection 
and transmission characteristics of plane waves at liquid-solid 
interface. The liquid is chosen to be inviscid and the solid 
half-space is homogeneous isotropic, thermally conducting 
viscoelastic. Both classical (coupled) and non-classical 
(generalized) theories of linear thermo-viscoelasticity have 
been employed to investigate the characteristics of reflected 
and transmitted waves. Reflection and transmission 

coefficients are obtained for quasi-longitudinal ( qP ) wave. 

The numerical computations of reflection and transmission 
coefficients are carried out for water-copper structure with 
the help of Gauss-elimination by using MATLAB software 
and the results have been presented graphically. 

Keywords- Reflection, Transmission, Viscoelastic Solid, 
Inviscid fluid, Critical angle. 

I. Introduction 

The problems of reflection and transmission of waves at 
an interface between liquid and solid media has many 
applications in under water acoustics and seismology. Ewing 
et al. [1], Hunter et al. [2] and Flugge [3] used mathematical 
models to accommodate the energy dissipation due to viscous 
effects in vibrating solids. Acharya and Mondal [4] 
investigated the propagation of Rayleigh surface waves in a 
Voigt-type [5] viscoelastic solid under the linear theory of 
non local elasticity. Schoenberg [6], Lockett [7], Cooper and 
Reiss [8] and Cooper [9] have investigated the problems of 
reflection and transmission of waves at an interface between 
viscoelastic isotropic media. 

In order to eliminate the paradox of infinite velocity of 
thermal signals in classical (coupled) thermoelasticity, Lord 
and Shulman [11] and Green and Lindsay [12] proposed 
nonclassical (generalized) theories of thermoelasticity which 
predict a finite speed for heat propagation. Sharma, etal. [13] 
studied the reflection of piezothermoelastic waves from the 
charge free and stress free boundary of transversely isotropic 
half space. 

In this paper, we discuss the reflection and transmission 
of plane waves at the interface between inviscid liquid half- 
space and thermo- viscoelastic solid half-space. The effects 
of incident angles and fluid loading on reflection and 
transmission coefficients are considered. The analytical 
results so obtained have been verified numerically and are 
illustrated graphically. 
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II. FORMULATION OF THE PROBLEM 

We consider a homogeneous isotropic, thermally 
conducting, viscoelastic solid in the undeformed state initially 

at uniform temperature T Q , underlying an inviscid liquid half 
space. We take the origin of the co-ordinate system 
(x, y, z) at any point on the plane surface (interface) and 
the z -axis pointing vertically downward into the solid half 
space which is thus represented by z > ■ We choose the x- 
axis along the direction of wave propagation in such a way 
that all the particles on the line parallel to the y-axis are equally 
displaced. Therefore, all the field quantities are independent 
of y -co-ordinate. Further, the disturbances are assumed to 
be confined to the neighborhood of the interface z—0 an d 
hence vanish as z — > °° ■ In the linear theory of 
homogeneous isotropic, the basic governing field equations 
of motion and heat conduction for solid and liquid (inviscid) 
medium, in the absence of heat sources and body forces, are 
given by 

p*V 2 u + (X + ft*) V V- U - J? V(T + tyS^ t) = pu ( 1 ) 
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Here A, jX are Lame's parameters, CC Q and a, arethermo- 

viscoelastic relaxation times and (X t is the coefficient of lin- 
ear thermal expansion, p is the density of the solid, 
T(x, Z, t) is the temperature change and 

u( X Z t) = (u w) ^ s tne displacement vector; fi[ is 
the thermal conductivity; C is the specific heat at constant 
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strain of the solid; t and t, are thermal relaxation times 
; /L is the bulk modulus, p and a * are the density and 
coefficient of volume thermal expansion, u is the velocity 
vector and T, is the temperature deviation in the liquid tem- 
perature from ambient temperature Tj,* ; 8 % is the 

Kronecker's delta with k — 1 forLS theory and fc — 2 for 
GL theory. 

The superposed dot notation is used for time 
differentiation. To facilitate the solution we define the 
following dimensionless quantities. 
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Here (q* is the characteristic frequency of the solid plate; 

£ is the thermomechanical coupling constant and c, , C 2 
are respectively, the longitudinal and shear wave velocities 
in the thermoelastic solid half-space; £ is the 

thermomechanical coupling and c, is the velocity of sound 

in the fluid. 

Upon using quantities (5) alongwith the relations 
u=#, x +ip, 2 , w = 4. K -1?, x (6) 

in equations (l)-(4), we get 
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in. BOUNDARY CONDITIONS 

The boundary conditions at z—0 can be expressed 

as T zz = -p , T xz = , w = w L , T_+H(T-T L ) = 
where H is the Biot's heat transfer constant. *• ' 

IV. SOLUTION OF THE PROBLEM 

We assume wave solutions of the form 



{0._ 5K, T. §}={A B._ C._D}ex.p{ik(xsm9-zcos8-cf)} 

(14) 

_ m 
where C — — is the non dimensional phase velocity, Q) is 
k 

the frequency and £ is the wave number. 

Upon using solution (14) in equations (8)-(12), we obtain 
a system of algebraic equations in unknowns A, B, C and D. 
The condition for the existence of non-trivial solution of this 
system of equations upon solving provide us 

k:=a: : co l t/ = l,2J>4) (15) 

where 

af + a\ = 1 - i&S£r 9 + ia? sr 9 z^l , 

fljfl, = -md a T a ,a 3 =— — -, a 4 =- I - - 

r = ^ +i&\ Tj = tfi* + to" 1 , 4 = t& + i a' 1 , 

al =cCq + ub>~ 1 , 

j5* = /? + i ty" 1 : of = a^ + 1 ty" 1 ;i 5 * = <5 + rfiT 1 

(16) 

In the absence of viscous effects ((Z =0 = 0.^ and 

thermal field (7" = = £ , 7" L = = £ L ) , we have 

1 1 

flj =L,fl 2 =r ,fl 3 =— a 4 =— (17) 

A. qP -Wave Incidence Upon a Plane Surface 

Let the suffix / and r represent incident and reflected waves, 

respectively. Omitting the term exp(— i (Ot) , we can write 

2 

= ft t -^ A rJ exp{ ikjixsm &. : + z cos 9j )} (1 S) 
W-Wt = -A.' h ^P 1 ! ^ (jc sin (9; + z cos ^ )} (19) 
T = T,=^S. : ,4 7 exp{fjt / (jcsin 9 j +zcos6^)} (2G) 

5^ — ^Li + ^r — J ^ ^P ^4 C-^ 5 ™ ^ + z COS0} + 

J 4 sxp{ik 4 (x^n0 4 —zcos& 4 )} (21) 
7^ =7j^ +T Lr =A exp{!'A" 4 (xsin i^+zcos0} + 

^^exp^Cxsin^-zcos^)} (22) 
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where 
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Upon using equations ( 1 8)-(22) in the boundary conditions 
(14) alongwith the fact that all the waves, incident, reflected 

and transmitted must be in phase at the interface g — for 
all values of x and t , we get 

jtj sin 8 l = k 2 sin 8 2 = k~ sin 8~ = k 4 sin 8 4 (24) 
This with the help of equation (15) implies that 

a x sin 8 X = a 2 sin 8 2 = fl? sin #, = o 4 sin 8 4 (25) 

The equation (25) is modified Snell's law in this situation. 
In the absence of thermal field, viscous effect and liquid, (25) 



becomes 



S sin 8 X = sin 8-. 



sin 8 X sin #, 
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The analytical expression of reflection and transmission 
coefficients /?* = A 4 I A t and 



■9^ _ 



A r / A ( (£ = 1,2,3) in the presence of thermal 



field for incident ^.P wave are obtained as 
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and A, , A 2 , A 3 , A. can be obtained from A. by replacing 
first, second , third and fourth column by 

[-a u a 24 a 34 0] respectively. 



here a n = COs2# 3 , 

a l3 = -sin2# 3 , 

a 2l = a*8 2 a 2 sin 2^, a 22 = a* 2 8 2 a 2 sm26 2 
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V. NUMERICAL RESULTS AND DISCUSSIONS 
In this section the reflection and transmission coefficients 
for qP wave incidence at an interface between thermo- 

viscoelastic solid and inviscid fluid have been computed 
numerically. The material chosen for this purpose is Copper, 
the physical data for which is given by Sharma, et al. [14] 

£ = 0.00265, 2 = 8.2xl0 10 AM- 2 , 



H = 4.2 xlO 10 Nm~ 



p = 8.950 xlO' kg m 



K = 1.13xl0 z Calm ~ L s l K ,a T =1.0x10 *K l , 
a =a x =6.883 1x10 13 ^, T =300K 

The liquid chosen for the purpose of numerical 
calculations is water, the velocity of sound in which is given 

by c L = 1.5x\0 3 m/s and density is 

p L =1000^m- 3 .r * = 29SK. 

Figs. 1 and 2 yields the behaviour of reflection / 
transmission coefficients for the angles of incidence of 
longitudinal wave propagates from fluid into solid. It is 
observed that for longitudinal wave incidence, the reflected 
longitudinal wave passes through a minimum at critical angles 

Q = 50° for elastic case which is known as Rayleigh-wave 

angle. At this angle a wave with large surface components is 
generated. These results parallel those obtained by Mott 
[10] in the analysis of incidence at a water-stainless steel 
interface, under the influence of dissipation. 
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Figure 1 . Qr -wave incidence at the interface (TVE/E) 
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Figure2. qP -wave incidence at the interface (VE/E) in the 
absence of thermal field 

CONCLUSIONS 

The reflection and transmission of plane waves at in viscid 
liquid- thermoviscoelastic solid interface has been analyzed 
theoretically. The significant effect of incident angle, thermal, 
viscosity and presence of liquid on the amplitude ratios of 
reflected and transmitted waves have been observed. 
Rayleigh angle phenomenon is explained. It is shown that 
reflected surface wave exist for incidence angles greater than 
Rayleigh wave angle as explained in [8]. 
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